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Thermodynamic, structural and electronic properties of isolated copper and iron atoms as well as their com- 
plexes with oxygen vacancies in tetragonal lead titanate are investigated by means of first principles calculations. 
Both dopants exhibit a strong chemical driving force for the formation of Mn — V& (M = Cu, Fe) defect asso- 
ciates. The most stable configurations corresponds to a local dipole aligned along the tetragonal axis parallel to 
the spontaneous polarization. Local spin moments are obtained and the calculated spin densities are discussed. 
The calculations provide a simple and consistent explanation for the experimental findings. The results are 
discussed in the context of models for degradation of ferroelectric materials. 
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I. INTRODUCTION 

Piezoelectric materials transform electric signals into me- 
chanical response and vice versa. They are widely used in 
sensors and actuators as well as in non- volatile memory de- 
vices. The properties of piezoelectric s and their behavior in 
the presence of electric fields are intimately linked to the de- 
fect structures (point defects, domain walls, grain boundaries) 
present in these materials. 

Defect dipoles, for instance, formed by the associa- 
tion of oxygen vacancies with metal impurities affect the 
switching behavior as well as aging and fatigue of piezo- 
electric materials J- Their existence has been shown both 
experimentally 1-4 and theoretically. 5 However, a comprehen- 
sive microscopic picture of the electronic, structural and ki- 
netic properties of these defect associates is still lacking. 
Such detailed knowledge is, however, a prerequisite for un- 
derstanding electrical "aging" and ferroelectric "fatigue" as 
exemplified by the model by Arlt and Neumann, 6,7 which 
relates these macroscopic degradation phenomena to micro- 
scopic processes. Since the time the first phenomenological 
models were developed^ quantum mechanical methods have 
become available which are both efficient and sufficiently reli- 
able for assessing the energetics of point defect and associates, 
directly. In the past, several studies addressed the properties 
of intrinsic point defects with particular attention to oxygen 
vacancies The role of dopants and impurities and the for- 
mation of defect dipoles has also been elucidated^r^ 4 - but im- 
portant aspects such as the dependence of the formation en- 
ergies on the Fermi level and the charge state have not been 
investigated in sufficient detail. Also, a comparative study of 
the different possible metal impurity-oxygen vacancy config- 
urations, which could serve as input for the defect models 
referenced above, is not available yet. Such information is, 
however, crucial both for understanding the generic properties 
of different impurities/dopants and for identifying properties 
which differentiate them, in particular with regard to their ef- 
fect on aging and fatigue. These considerations motivate the 
present work, in which density functional theory calculations 



are employed to study the energy surface for unbound oxygen 
vacancies and oxygen vacancies complexed with Fe or Cu im- 
purities in tetragonal lead titanate. The results are discussed 
in the context of existing defect models and have been used to 
interpret recent electron spin resonance experiments. 15 



II. METHODOLOGY 
A. Computational details 

Quantum mechanical calculations based on density func- 
tional theory (DFT) were carried out using the Vienna 
ab-initio simulation package (VASP) which implements the 
pseudopotential-plane wave approach i^"— The projector- 
augmented wave method was employed for representing the 
ionic cores and core electrons j 2 ^- In accord with earlier cal- 
culations (see e.g., Refs. [8| andHj) the bd electrons of Pb and 
the 3 s and 3p electrons of Ti were treated as semi-core states, 
while for Cu and Fe the valence configuration included the 4d, 
3d and 3p electrons. All calculations were performed using 
Gaussian smearing with a width of a = 0.2 eV and a 4 x 4 x 4 
Monkhorst-Pack mesh for Brillouin zone sampling. 22 Calcu- 
lations were carried out using the local spin density approxi- 
mation (LSDA) to account for the magnetic behavior of iron. 50 
Some calculations were also repeated using the non-spin po- 
larized local density approximation (LDA). The difference be- 
tween the LDA and LSDA results allows to estimate the spin 
contribution for a particular defect or charge state. 

With this computational setup, we obtain a lattice constant 
of do = 3.866 A and an axial ratio of c/a = 1.05 for the 
tetragonal phase of PbTiC>3 in reasonable agreement with ex- 
periment. The calculated band gap of 1 .47 eV is considerably 
smaller than the experimental value but consistent with the 
well known LSDA band gap error. Since there is no unique 
procedure for correcting the band gap error (see e.g, Refs.l23l- 
\25b . in the present work, we abstain from any explicit cor- 
rections and discuss the electronic transition levels relative to 
the calculated band edges instead (see also our discussion in 
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FIG. 1. (Color online) Finite-size scaling of formation enthalpies 
to obtain formation energies in the limit of infinite dilution for three 
of the most important defects. The small filled and large open cir- 
cles show the data before and after application of the monopole- 
monopole correction term. 



secima). 

In order to model the defect configurations, super cells con- 
taining 2x2x2 to 2x2x4 unit cells equivalent to 40 to 80 
atoms were employed. Each configuration was relaxed until 
the maximum force was less than 15 meV/A. For charged de- 
fects a homogeneous background charge was added. Finite- 
size scaling was employed in order to obtain the formation 
and binding energies at infinite dilution i2Si2Z For the charged 
defect calculations monopole-monopole interactions were ex- 
plicitly corrected 2 ^ 2 ^ using values of ef x = 100 and ef 3 = 34 
for the static dielectric tensor, 51 while higher order terms were 
implicitly included by the finite- size scaling procedure j2&2Z 
Both the effect of the monopole-monopole term and the finite- 
size procedure are illustrated in Fig.[T]for three of the most im- 
portant defects which demonstrates that the finite- size scaling 
approach works reliably in the present case. The extrapolation 
errors of the finite- size scaling procedure are given in Tables H 
HH andim These values are a measure for the quality of the 
finite- size scaling procedure only, and do not include other er- 
rors intrinsic to DFT. 



B. Formation energies 

The formation energy of an intrinsic defect in charge state, 
q, depends on the relative chemical potentials of the con- 
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where Ed is the total energy of the defective system, Eh is 
the total energy of the perfect reference cell, £Vbm is the po- 
sition of the valence band maximum, Arii denotes the differ- 
ence between the number of atoms of type i in the reference 
cell with respect to the defective cell, and ii\ ulk is the chemical 
potential of the reference phase of atom type i. In the present 
context we are not concerned with the effect of the chemi- 
cal environment (Afii in equation (Q])) but focus on the transi- 
tion levels and the binding energies between oxygen vacancies 
and metal impurities/dopants (see equation Q below), both of 
which are independent of the chemical potentials. 

The binding energy for the metal impurity-oxygen vacancy 
(Mri — Vq) associates considered below is calculated as the 
difference between the formation energies of the defect com- 
plex and the isolated defects 

E b = AE f D [M Tl - V ] - AE f D [M Ti ] - AE D [V ] . (2) 

The binding energy depends on the Fermi level due to the 
Fermi level dependence of the formation energies. For a given 
value of ii e the difference in equation © is calculated between 
the most stable charge states for each defect involved. By in- 
sertion of equation (Q]) into © one can show that the terms 
with Arii cancel, whence E b is independent of the chemi- 
cal potentials, A/^. Following the sign convention introduced 
via equation ([2]), negative values for E b imply that the system 
gains energy by association, whereas positive values imply an 
effective repulsion between the isolated defects. 



III. RESULTS 

A. Isolated oxygen vacancies 

The crystallographic structure of a tetragonal perovskite 
possesses two symmetrically inequivalent oxygen sites which 



TABLE I. Formation energies (eV) for unbound oxygen vacancies 
in tetragonal lead titanate for a Fermi level at the valence band max- 
imum (/x e = OeV) and metal-rich conditions (A/ipb = Ajin = 
eV). The finite- size scaling extrapolation errors are given in brack- 
ets in the last column. 
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FIG. 2. (Color online) Variation of the formation energy of unbound 
oxygen vacancies with the Fermi level for extreme metal-rich con- 
ditions (A/xpb = A/ixi = OeV). The white area corresponds to the 
calculated band gap. The numbers indicate the defect charge state. 



are located within the afr-plane and along the c-axis, respec- 
tively. Based on DFT calculations, Park and Chadi 8 identi- 
fied one possible vacancy configuration along the c-axis (V c ) 
and two distinct vacancies configurations within the afr-plane 
(V*™, V^ b d ). The latter two configurations differ in the orien- 
tation of the polarization above and below the vacancy plane. 
In the present work we have found the "up-down" configura- 
tion (V^jf, see Ref. [8| for nomenclature), in which the polariza- 
tion vectors are oriented head-to-head in the vacancy plane, to 
be highly unstable with respect to the "switchable" configura- 
tion (V^™), which is characterized by a uniform orientation of 
the polarization vector. Although Park and Chadi did not pro- 
vide a quantitative energy comparison for and V^jf due to 
missing data for the 180°-domain wall energy, they estimated 
the V*™ configuration to be lower in energy.- For these rea- 
sons, we focus on the V£™ configuration which for the sake 
of brevity in the following is denoted V a b> 

Table U summarizes the formation energies for unbound 
oxygen vacancies obtained for extreme metal-rich conditions 
(A/ipb = A/iTi = eV). In all charge states the c-type va- 
cancy is energetically preferred over the afr-type vacancy in 
agreement with the calculations in Ref. 8. For infinite dilution 
(extrapolated values in TableUthe energy difference is as large 
as 0.8 eV for q = +2 which is on the order of magnitude of 
the energy barrier for migration^ 

The variation of the formation energy with the Fermi level 
is shown in Fig. [2] Practically over the entire band gap oxygen 
vacancies are doubly positively charged. For the given chem- 
ical conditions the formation energies are negative which in- 
dicates that the material is unstable for these chemical poten- 
tials (A/ipb = A/ixi = eV). If the oxygen partial pressure is 
raised (i.e., A/iTi and A/i P b become less negative while A/i0 
becomes more negative), the oxygen vacancy formation en- 
ergy increases and the stability condition is fulfilled. 

On first sight the present finding is at variance with the 
calculations in Ref. 1 1 in which the +2/0 transition level of 
the oxygen vacancy is located in the middle of the band gap. 



FIG. 3. (Color online) Fermi level dependence of formation en- 
ergies of free and complexed Cu impurities. Parallel lines imply 
equal charge states as indicated in the plot. The white stripe indi- 
cates the (calculated) band gap while the gray stripes indicate the 
valence and conduction bands, respectively. Extrapolated formation 
energies from Table HTl were used except for Cun — vj 4 ^ for which 
the data obtained with 80-atom cells are shown. 

This finding is, however, a consequence of the corrections ap- 
plied in that work: In order to correct for the the band gap 
error the conduction band minimum was rigidly shifted up- 
wards without correcting the formation energies. If this shift 
is omitted the results of Ref. fill and the findings of the present 
work are consistent with each other. The approach in Ref. [U 
ignores the fact that in the case of singly charged (V^) and 
neutral oxygen vacancies (V^ ) the excess electron(s) occupy 
conduction band states. In order to maintain internal consis- 
tency, a correction term 2 - on the order of z e AEc would have 
to be added to the formation energies of Vq and V@ , where z e 
is the number of occupied conduction band states and AEq 
denotes the difference between the calculated and the exper- 
imental band gap. 52 If these energy terms are included, the 
+2/+1 (+2/0) transition is no longer located in the middle 
of the (experimental) band gap but near the conduction band 
minimum. This is equivalent to the location of the transition 
level with respect to the calculated CBM if no correction is 
applied. 

1. Isolated copper impurities 

Since there is experimental evidence that copper impuri- 
ties in lead titanate preferentially occupy Ti-sites^ 2 - - — only the 
substitutional position was considered (Cun) in the present 
calculations. The calculated formation energies for different 
charge states of CuTi are compiled in Table HB for extreme 
metal-rich conditions (A/ipb = A/iTi = A/ic u = A/i Fe = 
eV). They can be used to derive the Fermi level dependence 
of the defect formation energy as demonstrated in Fig. [3] which 
shows that uncomplexed Cu-impurities can occur in charge 
states q = — 1 and —2. The transition between these two 
charge states is located ~ 0.2 eV below the middle of the cal- 
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TABLE II. Formation energies (LSDA) in units of eV for isolated 
and complexed Cu-impurities on Ti- sites (Cuxi) in tetragonal lead 
titanate for a Fermi level at the valence band maximum (/x e = eV) 
and metal-rich conditions (A/ipb = A /in = A/icu = OeV). The 
finite- size scaling extrapolation errors are given in brackets in the last 
but one column. An] [ denotes the difference between the number of 
electrons in spin-up and spin-down states. 
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culated band gap. 

In the defect chemistry of oxides, often a fully ionic pic- 
ture is adopted e.g., in Cu-doped lead titanate Cu 2+ ions are 
usually assumed to replace Ti 4+ ions with two defect elec- 
trons being associated with the defect site. In Kroger- Vink 
notation the resulting defect is written as Cu Ti . Accordingly, 
copper-impurities should act as acceptors and create holes in 
the valence band. 

The present quantum-mechanical calculations are free of 
a-priori assumptions with respect to the oxidation state of a 
defect: Different charge states, q, are obtained by varying the 
number of electrons in the entire system. In order to provide 
a connection between the ionic picture outlined above and the 
charge states given in Table [II] and Fig. [5J one must therefore 
analyze the electronic structure of the defect. 

The partial and total densities of states for a neutral (q = 0) 
copper impurity in comparison with the ideal (defect-free) 
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FIG. 4. (Color online) Spin-decomposed total and partial density 
of states for neutral uncomplexed Cu impurities in lead titanate. The 
partial densities of states for the Cu atom have been scaled by a factor 
of twenty and are shown separately in the upper part of the figure. 
The density of states for the ideal crystal is shown below. In the 
latter case the density of states is independent of the spin orientation 
A Gaussian filter with a — 0.05 eV was applied to optimize the 
visualization. 



system are shown in Fig. |4l The bottom of the conduction 
band is predominantly composed of Ti-3d orbitals (feature A 
in Fig. [4]). Since in this energy range the contribution Cu-3d 
orbitals is practically zero, copper substitution leads to an ef- 
fective reduction of the density of states at the bottom of the 
conduction band. With regard to the top of the valence band, 
the most obvious feature is a hybridization of Cu-3d orbitals 
with 0-2p orbitals (feature B in Fig. [4]) while the contribu- 
tion of Cu-4s and Cu-3p states is negligible. In addition, a 
Cu-induced state emerges attached to the valence band edge 
(feature C in Fig.@]). As electrons are added to the system, this 
defect induced state is, however, shifted into the band gap be- 
cause the excess charge is localized in a small volume around 
the defect. As a result the local potential at the copper atom 
is modified and an energy offset of the electronic states of this 
defect is induced. 

Integrating the density of states and the occupied levels up 
to the valence band edge and including the defect induced 
peak, reveals three empty levels (holes) for the neutral charge 
state. As electrons are added to the system (q < 0) these 
empty states are gradually filled. For the charge state — 1 there 
are two holes localized at the Cu atom. The situation is thus 
equivalent to placing a Cu ++ on a Ti-site giving Cu Ti , while 
for q = — 2 there is one unoccupied valence band level corre- 
sponding to Cu Ti (Cu + on Ti 4+ ). If the Fermi level is located 
in the middle of the band gap, the most stable defect is Cu Ti . 

The — 1 / — 2 transition can be understood as a local reduc- 
tion of copper: The transition from charge state q = —1 to 
q = —2 is equivalent to the reaction Cu ++ + e~ —> Cu + . 
Experimentally, the redox potential for this reaction has been 
measured in aqueous solution to be +0.16 eV with respect to 
the standard hydrogen electrode, which on an absolute energy 
scale is located at 4.5 eVM Thus, the Cu ++ /Cu + redox po- 
tential is 4.7 eV. On the other hand, the electron affinity of 
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FIG. 5. (Color online) Isolated copper impurity in charge state 
q = — 1 (a-d) and copper impurity-oxygen vacancy complex in 
charge state q = (e-h). (a,e) Position of atoms in the (100)- 
plane containing the defect, numbers given in square brackets give 
the equivalent distances in the defect-free crystal; (b,f) total charge 
density, a logarithmic scale has been chosen for the charge density in 
order to enhance the features, (c,g) difference between total charge 
density of defective cell and defect-free cell, and (d,h) spin densities. 
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FIG. 6. (Color online) Schematic of different conceivable vacancy- 
impurity atom configurations. The tetragonal perovskite lattice is 
shown along the [100] axis and Pb atoms have been omitted for clar- 
ity. The different gray shaded areas indicate the size of the different 
supercells employed in the present work. 



ison with the ideal crystal, copper substitution causes the B- 
site to move towards the oxygen (100)-plane. For the charge 
state q = — 1 (q = — 2) the relative displacement of the copper 
atom with respect to its in-plane oxygen neighbors is 0.15 A 
(0.19 A) along [001], compared to 0.28 A in the defect free 
crystal. Copper substitution thus causes a relaxation of the 
B- site atom in the direction of the nearest oxygen plane. The 
charge density plots show that the defect charge (Fig. [3c)) as 
well as the spin-asymmetry (Fig. (3d), also compare Fig. (4j> 
is localized at the copper atom. In addition, the difference 
between the spin-densities in the (100) -plane reveals a slight 
asymmetry with respect to the tetragonal axis. Analysis of the 
electron and spin density at farther atoms - in particular the 
nearest Pb atoms in the (110) -plane - indicates that both of 
them are virtually unaffected by the presence of the Cu impu- 
rity. 



PbTi0 3 is 3.5 eV and the band gap is about 3.4 eV (Ref.[H 
and references therein) placing the valence band maximum 
(VBM) 6.9 eV below the vacuum level. In analogy, the re- 
dox reaction occurs roughly in the middle of the band gap 
(1.2 eV below the CBM or 2.2 eV above the VBM), which is 
in good agreement with the present calculations which locate 
the — 1/ — 2 (or Cu ++ /Cu + ) transition level near the mid gap. 
This analogy has, however, merely a qualitative character. 

In order to obtain a measure for the spin-polarization of 
copper defects we calculated the difference between the num- 
ber of electrons in spin-up and spin-down states (An^ [ in 
Table [II|). In both relevant charge states (q = — 1 and —2) 
unpaired electrons are present in the system which allows to 
probe the defect center spectroscopically. 

The geometry and charge density of an isolated copper im- 
purity in charge state q = — 1 is shown in Fig.[5j In compar- 



2. Copper-vacancy complexes 

As illustrated in Fig. [6| different configurations for defect 
associates of oxygen vacancies with 5-site impurities are con- 
ceivable. Since the impurity atom breaks the translational 
symmetry along the c-axis, the degeneracy of the positions 
above and below the impurity atom is lifted. In the first 
neighbor shell one can distinguish configurations in which 
the oxygen vacancy is placed along the c-axis either below 
(Mri — Vc X \ shorter Mn — distance) or above the im- 

(2) 

purity atom (M T i — Vc , longer M T i — distance), or on 
one of the four symmetry equivalent sites within the afr-plane 
(Mri — V a b)> If one considers the second-nearest neighbor 
shell, another two configurations involving c-type vacancies 
are possible (Mri — Vq 3 \ Mri — V^). Further configurations 
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FIG. 7. (Color online) Binding energies for Cu impurities in tetrag- 
onal lead titanate obtained from the extrapolated formation energies 
given in Tables U and HH Negative energies imply a chemical driving 
force for association. 



are obtained by placing the vacancy on second-nearest neigh- 
bor afr-sites (e.g., M Ti - V$ and M Ti - V$ in Fig.©. In 
the present work, we have considered first and second-nearest 
neighbor configurations involving c-type vacancies as well as 
first-nearest neighbor configurations with ab- vacancies. Ex- 
ploratory calculations were also carried out for second-nearest 
neighbor afr-type vacancy associates. This possibility was, 
however, not pursued further, since the total energies were 
significantly larger than for the respective lowest energy con- 
figuration. 

The Curi — Vc 1 ^ complex, in which the vacancy is located 
at the nearest oxygen site along the c-axis (compare Fig. (6[ 
and Fig. 0), is the most stable configuration. The most sta- 
ble charge state is q = as expected if one combines V" (see 

Fig. [2]) and Cu Ti . The energy difference with respect to the less 

stable complexes ([Curi — V a b] x and [Curi — V^] x ) is rather 
large (0.8-1.0 eV). This large energy difference is particu- 
lar noteworthy since previous model calculations assumed a 
much smaller energy difference on the order of 0.06 eV (com- 
pare Ref.[Z|). 

It is very instructive to compare the charge densities for iso- 
lated and complexed copper impurities. As shown in Fig. [3d), 
the spin density for the isolated copper impurity is symmet- 
ric with respect to four-fold rotations about the [001] axis. 
This symmetry also applies for the copper- vacancy complex 
[Fig. Oh)] but in addition the spin density pattern contains a 
(001) mirror plane through the Cu site. As can be seen by 
comparison of Figs. [3a) and (e) this distinction results from 
the shift of the Cu atom which in the presence of a vacancy 
on the nearest O-site relaxes almost completely into the (001) 
plane of oxygen atoms leading to a local pseudo-cubic en- 
vironment. These observations have important implications 
for the interpretation of electron paramagnetic resonance mea- 
surements as discussed in Ref . [l5|- 

Finally, one can use the formation energies of isolated oxy- 
gen vacancies, copper impurities and their complexes to de- 



FIG. 8. (Color online) Fermi level dependence of formation en- 
ergies of free and complexed Fe impurities. Parallel lines imply 
equal charge states as indicated in the plot. The white stripe indi- 
cates the (calculated) band gap while the gray stripes indicate the 
valence and conduction bands, respectively. Extrapolated formation 
energies from Table Hill were used except for Fexi — Vc^ for which 
the data obtained with 80-atom cells are shown. 

rive the binding energy and its dependence on the Fermi level. 
As shown in Fig. [7J there is a very strong chemical driving 
force for association of copper impurities with oxygen vacan- 
cies irrespective of the Fermi levels. If there are more oxygen 
vacancies than copper impurities in the system, all copper im- 
purities will be complexed. In the opposite scenario all oxy- 
gen vacancies would be complexed and the number of isolated 
copper impurities would be reduced by the number of metal 
impurity-oxygen vacancy associates. 



3. Iron impurities 

The formation energies for uncomplexed Fe impurities are 
compiled in Table [III1 and presented as a function of the Fermi 
level in Fig. [8l The uncomplexed iron impurity occurs pre- 
dominantly in the neutral charge state (q = 0). Transition lev- 
els are only present near the band edges. It is noteworthy that 
according to Fig. [8] uncomplexed Fe should display ambipo- 
lar behavior. Analysis of the total and partial density of states 
shows that iron induces a defect level which for the neutral 
charge state is located in the middle of the band gap (feature 
A in Fig. [9]). The density of states for the two different spin 
orientations differ considerably. While one spin-orientation 
(spin-down in Fig. [9]) gives rise to a level in the band gap (fea- 
ture B in Fig. [9]), the other one leads to additional states at 
the top of the valence band (feature C in Fig. 0. As in the 
case of copper, strong hybridization occurs between the im- 
purity atom 3d-orbitals and the 0-2p levels of the host mate- 
rial. Also, in equivalence to copper Fe-3d levels do practically 
not contribute to the conduction band, slightly diminishing the 
density of states at the bottom of this band. Due to the pres- 
ence of gap states defect induced holes and electrons cannot 
be unambiguously identified. In contrast to copper, a simple 
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TABLE III. Formation energies (LSDA) in units of eV for isolated 
and complexed Fe-impurities on Ti- sites (Fen) in tetragonal lead ti- 
tanate for a Fermi level at the valence band maximum (/i e = eV) 
and metal-rich conditions (A/xpb = A /in = A/iFe = OeV). The 
finite- size scaling extrapolation errors are given in brackets in the last 
but one column. Art] [ denotes the difference between the number of 
electrons in spin-up and spin-down states. 



Defect 


An 1t 


40 


60 


80 


extrapolated 


Fen 














-2 





8.64 


8.81 


8.87 


9.11 


(0.04) 


-1 


1 


7.14 


7.18 


7.21 


7.27 


(< 0.01) 





2 


6.26 


6.23 


6.23 


6.19 


(0.01) 


+ 1 


3 


6.15 


6.12 


6.10 


6.05 (< 0.01) 


Fen - 


V ab 












-2 





8.69 


8.91 


8.96 


9.25 


(0.07) 


-1 


1 


6.94 


6.99 


7.01 


7.08 (< 0.01) 





2 


5.46 


5.44 


5.44 


5.41 


(0.02) 


+ 1 


3 


4.49 


4.43 


4.40 


4.31 (< 0.01) 


+2 


4 


4.18 


4.14 


4.09 


4.01 


(0.04) 


Fen - 














-2 





8.42 


8.52 


8.53 


O CI K 
O.00 


(U.UoJ 


-1 


1 


6.50 


6.52 


6.53 


6.56 


(< 0.01) 





2 


5.00 


4.93 


4.93 


4.84 


(0.04) 


+ 1 


3 


3.99 


3.88 


3.85 


3.70 


(0.02) 


+2 


4 


3.82 


3.73 


3.68 


3.54 


(0.02) 


Fen — 














-2 





8.61 


9.08 


9.43 


10.20 


(0.13) 


-1 


1 


6.83 


7.17 


7.39 


7.92 


(0.06) 







5.40 


5.60 


5.46 


5.62 


(0.27) 


+ 1 




4.34 


4.25 


4.25 


4.15 


(0.05) 


+2 


4 


4.09 


4.16 


4.17 


4.25 


(0.03) 


Fen — 


K (3) 












-2 







9.19 


9.11 


8.86 




-1 







6.97 


6.87 


6.55 












5.26 


5.21 


5.07 




+ 1 


1 




4.32 


4.23 


3.96 




+2 


2 




4.24 


4.10 


3.69 




Fen - 


K (4) 












-2 









9.42 






-1 









7.21 
















5.55 






+ 1 


1 






4.54 






+2 


2 






4.27 







a The difference between electrons in spin-up and spin- 
down states, Ari] [, for the configuration with the lowest 
energy is system size dependent. 

correlation between the Kroger- Vink notation for defects and 
the charge states, q, is, therefore, not possible. 

The potentials for the redox reactions Fe 3+ + e~ — >• Fe ++ 
and Fe ++ + 2e~ —> Fe° have been measured in aqueous so- 
lution to be 0.77 eV and — 0.41 eV, respectively. Based on 
these values and using the data for the electron affinity and 
the band gap cited above, one would expect transition levels 
at 1.6 eV (Fe 3+ /Fe ++ , i.e. Fq^/Fq^) and 2.8 eV (Fe ++ /Fe°, 
i.e. Fe Ti /Fe^) with respect to the experimental valence band 




-8 -6 -4 -2 2 

Energy (eV) 



FIG. 9. (Color online) Spin-decomposed total and partial density 
of states for neutral uncomplexed Fe impurities in lead titanate. The 
partial densities of states for the Fe-d orbitals are shown separately 
in the upper part of the figure. The Fe-s and Fe-p states have been 
omitted since in the energy window shown they are practically zero. 
A Gaussian filter with a = 0.05 eV was employed to smoothen the 
data. 



maximum. In contrast, the present calculations indicate tran- 
sitions about 0.15 eV above the VBM and 0.39 eV below the 
CBM. The analogy with the behavior of ions in solution, 
which worked reasonably in the case of copper, thus appears 
to fail for iron. 



4. Iron-vacancy complexes 

The iron-vacancy complexes which are analogous to the 
copper- vacancy associates described in Sec. lIII A2l and Fig.[6l 
In agreement with earlier calculations^ the configuration with 
the shortest distance between impurity center and oxygen va- 
cancy (Fexi — Vq 1 ^) is the most stable (Fig. [8]). The energetic 
ordering of the remaining configurations is, however, some- 
what different from the case of copper. This observation has 
implications for the migration of oxygen vacancies and is dis- 
cussed in the context of point defect models in Sec. [TV] and 
Ref.HH 

The binding energies for Fe impurities shown in Fig. [TT] 
are in general smaller than for copper and for Fermi levels 
very close to the valence band the binding energy can even 
become positive. Nonetheless for most conditions there is a 
rather strong (^> ksT) driving force for association. 

Analysis of the spin density shows that for both isolated 
and associated iron impurities unpaired electrons are local- 
ized in the d-electron states at the iron site (Table [Till). It fur- 
thermore reveals that in both cases the spin density pattern is 
symmetric with respect to fourfold rotations about the tetrag- 
onal axis. Unlike in the case of copper- vacancy associates the 
(001) plane is, however, not a mirror plane. 
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FIG. 10. (Color online) Isolated iron impurity in charge state q = 
(a-d) and iron impurity-oxygen vacancy complex in charge state 
q = +1 (e-h). (a,e) Position of atoms in the (lOO)-plane contain- 
ing the defect, numbers given in square brackets give the equivalent 
distances in the defect-free crystal; (b,f) total charge density, a loga- 
rithmic scale has been chosen for the charge density in order to en- 
hance the features, (c,g) difference between total charge density of 
defective cell and defect-free cell, and (d,h) spin densities. 



IV. DISCUSSION 

The results obtained in this study need to be discussed in the 
context of degradation phenomena in ferroelectric materials. 
Two distinct processes have been widely described: (1) The 
gradual degradation of ferroelectric properties (usually over 
extended periods of time) in the absence of electric fields is 
termed "aging",— This process is commonly attributed to an 
increasing restriction of domain- wall motion with time— and 
typically accompanied by a shift of the ferroelectric hystere- 
sis loop along the electric field axis ^242 The shift is known as 
the internal bias field and it has been proposed to be related to 
the relaxation of defect dipoles £i In general, acceptor-doped 
ferroelectric ceramics with high concentrations of oxygen va- 
cancies are particularly prone to aging. (2) The degradation 
of ferroelectric properties in the presence of a (usually oscil- 
lating) electric field is termed electric "fatigue". There are 
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FIG. 11. (Color online) Binding energies as defined by equation 
© for Fe impurities in tetragonal lead titanate obtained from the 
extrapolated formation energies given in Tables U and |ni| Negative 
energies imply a chemical driving force for association. 



various contributions to this phenomenon and possible expla- 
nations involve defects of all kinds. 

The most widely known model for aging 6,7 relates the ap- 
pearance of internal bias fields to the existence and gradual re- 
orientation of impurity-oxygen vacancy defect dipoles. Sim- 
ilar models have been discussed in the literatur e^ 41 ! 42 Exper- 
imentally, an increased internal bias field has been observed 
for aged lead zirconate-titanate compounds and was attributed 
to defect dipole alignment leading to a clamping of domain 
walls. 43 

The availability of microscopic information on the energy 
landscape for defect formation and migration is the key ingre- 
dient for the these models. The present work provides a data 
base for these model and allow to verify their basic assump- 
tions. 

The key information regarding the orientation and the mag- 
nitude of the defect dipoles as well as their energetic ordering 
is summarized in Fig. [12] This picture needs to be compared 
to the energy landscape proposed by Arlt and Neumann^ on 
the basis of electrostatic arguments (Figs. 4 in Refs. @ and 
0). The latter analysis was carried out for barium titanate and 
involved several crucial assumptions with respect to e.g., the 
magnitude of the local dipole moment or the dielectric con- 
stant in the vicinity of the defect center. In agreement with 
the conclusions of the Arlt-Neumann model, the quantum- 
mechanical calculations show that for the energetically most 
stable configuration (Mn — Vc 1 ^) the defect induced dipole 
is oriented parallel to the overall spontaneous polarization, P s 
(compare Fig. [T21>. Thus in equilibrium (i.e., a state which 
would be achieved after "infinite" aging) the polarization lo- 
cally increases due the formation and alignment of local de- 
fect dipoles. In contrast, the defect dipole of the alternative 

(2) 

c-axis configuration (Mn — Vc ) causes a local reduction of 
the polarization. 

In Ref . two estimates for the energy differences between 
different dipole orientations are derived. The first one treats 
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FIG. 12. (Color online) Summary of key results for the most important charge states of nearest neighbor Cu/Fe- vacancy complexes ([Cun — 
Vjz>] x , [Fexi — Vja]"). The large arrow on the far left indicates the direction of the spontaneous polarization with respect to the ideal unit cell. 
The short arrows indicate the direction of the defect induced dipole moments. Distances are given between the impurity atom and the position 
of the vacancy (as defined by the positions of the neighboring oxygen planes. Values in square brackets refer to the ideal (undistorted) lattice. 
The lower panel visualizes the formation energy differences. 



the defect as a sphere with a defect induced polarization em- 
bedded in a ferroelectric matrix. It yields values on the or- 
der of 0.5 to 1.2 eV (compare equations (8a-8c) in Ref. 0) 
which are of a similar magnitude as the values obtained from 
the present calculations. The second estimate is based on the 
model of an ideal point dipole and yields much smaller en- 
ergy differences on the order of 0.03 to 0.06 eV (equation (9) 
in Ref. 7). It is important to realize that the modeling of aging 
carried out in the same paper was actually based on the latter 
values. 



In fact the energy differences obtained in the present work 
are actually on the order of magnitude of the barriers for the 
migration of free oxygen vacancies in this system (0.8-1.0eV, 
see e.g. Ref.lTTl). In the vicinity of impurities the energy sur- 
face, which describes the barriers for the transformation of 
different defect dipole configurations into each other, is there- 
fore severely distorted^ This effect is important for the be- 
havior of defect dipoles in the presence of oscillating electric 
fields, which pertains to the fatigue of the material. 



Finally, Fig . [12] highlights a noteworthy difference between 
Cu and Fe. While for Cu-complexes the energy difference 
between [Curi — V^] x and [Cun — V^] x is as large as 
1.21 eV, for [Fe Ti - V c W ] m and [Fe Ti - V c {2) ] 9 the difference 
is just 0.45 eV. This difference has implications for the redis- 
tribution of the defect dipoles in the presence of an oscillating 
electric field and thus for the process of fatigue. A detailed 
account of these behavior will be published elsewhere 



V. CONCLUSIONS 

This paper presents a comprehensive and detailed investi- 
gation of Fe and Cu dopants/impurities in tetragonal lead ti- 
tanate. Both metal ions are found to be very effective traps 
for oxygen vacancies. In the most stable configurations the 
oxygen vacancy is located on a c-site such that the Mn — V® 
distance is minimal, and the defect dipole is aligned paral- 
lel the spontaneous polarization. This result is in qualitative 
agreement with an electrostatic analysis of the defect induced 

(2) 

dipole field. An alternative configuration (M T i — Vc ) in 
which the dipole is oriented anti-parallel with respect to the 
spontaneous polarization is somewhat higher energy but is an- 
ticipated to be of importance for the fatigue of lead titanate 
based ferroelectrics. 

Since all relevant defects considered in this investigation 
carry unpaired electrons, they are visible for electron spec- 
troscopic methods. In fact, by means of analyzing the spin 
density patterns for different defects it was possible to in- 
terpret recent electron paramagnetic resonance measurements 
and to resolve an apparent disagreement between experiment 
and calculations.-^ The results obtained in this study represent 
important information on the energy landscape for energy for- 
mation and association and provide the basis for the improve- 
ment of defect models for the degradation of ferroelectric ma- 
terials. 
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